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Abstract 

It appears well accepted in the literature that the correlator of Polyakov 
loops in a finite temperature system decays with the "average" free energy of 
the static quark-antiquark system, and can be decomposed into singlet and ad- 
joint (or octet for QCD) contributions. By fixing a gauge respecting the transfer 
matrix, attempts have been made to extract those contributions separately. In 
this paper we point out that the "average" and "adjoint" channels of Polyakov 
loop correlators are misconceptions. We show analytically that all channels re- 
ceive contributions from singlet states only, and give a corrected definition of the 
singlet free energy. We verify this finding by simulations of the 3d SU(2) pure 
gauge theory in the zero temperature limit, which allows to cleanly extract the 
ground state exponents and the non-trivial matrix elements. The latter account 
for the difference between the channels observed in previous simulations. 



1 Introduction 



In recent years there has been a growing interest in assessing possible contributions of 
color octet bound states to physical processes. While such states carry color charge and 
hence do not exist as asymptotic particle states, they may appear as virtual intermedi- 
ate states in certain processes to whose cross sections and branching ratios they should 
then measurably contribute. Examples are quarkonium production where such contri- 
butions arise naturally in the framework of Non-Relativistic QCD pp, and the QCD 
plasma phase in which colored states might be excited thermally 0. It is expedient to 
first try and understand such physics in the heavy quark limit, where bound states can 
be obtained from a non-relativistic Schrodinger equation. This requires knowledge of 
the non-perturbative static quark- ant iquark potential as input, which may be extracted 
from lattice simulations of Wilson or Polyakov loops. 

While the Wilson loop probes the color singlet potential, there seems to be a common 
understanding in the literature that the correlator of Polyakov loops decays exponen- 
tially with the "color average" potential, which decomposes into two parts: one related 
to the color-singlet (1) potential and one to the color-octet or adjoint (A) potential 
0111 IS]- These potentials are thought to exist, and be distinct in the zero-temperature 
limit, corresponding to different sectors in the Hilbert space with color sources coupled 
to the singlet and adjoint representations, respectively. This is corroborated by per- 
turbation theory where the zero-temperature potentials can be computed to low orders 
in the coupling and are found to be gauge-invariant. To leading order one has 015] 

(N 2 -l)V A {r) = -V l {r) + 0{g i ). (1) 

More recently it was shown that, by dressing (or gauge fixing) the Polyakov loop in 
a way that respects the transfer matrix, it is possible to extract the singlet potential 
Vi(r) from Polyakov loop correlators jB], which suggests to obtain the adjoint channel 
by subtracting the singlet piece from the average. This approach was used in several 
numerical investigations of the different color channels at finite temperature [7], for 
which the potentials turn into superpositions of Boltzmann weighted excitations and 
are interpreted as free energies. A recent simulation employing different gauges con- 
firmed gauge independence for the zero temperature potentials, but observed gauge 
dependence at finite temperatures jH]. 

In this paper we clarify the situation regarding Polyakov loops. Converting Euclidean 
expectation values into traces over states in the Hilbert space, we show analytically 
that Polyakov loop correlators in all channels receive contributions from singlet states 
exclusively. Thus none of the channels or their combinations can serve to define an octet 
potential. Moreover, at non-zero temperature the standard Polyakov loop correlator 
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does not average over the color channels, but probes a thermal sum of only singlet 
excitations. 

We verify these statements by numerical simulations of the pure SU(2) gauge theory 
in 2+1 dimensions in the zero temperature limit. We compute the correlators in all 
channels at a series of low temperatures and extract the energy of the lowest state as 
well as its weight in the thermodynamic sum (its matrix element). The energies are 
identical within errors for all channels. We further show that the putative difference 
between the color channels, and in particular the repulsive r-dependence of the "ad- 
joint" potential observed in previous simulations, is actually due to the r-dependence 
of the matrix element of the singlet ground state. 



2 Operators for the different channels 

On an L 3 x N t lattice with periodic boundary conditions, we consider the untraced 
Polyakov loop (timeslices are labeled from O...N t — 1) 

Nt-l 

L(x)= J] tfo(x,t). (2) 

t=0 

According to jSllHEl, the color average potential at temperature T = l/(aN t ) is defined 

by 

-VMr,T)/T = _l (Tr L t (x)TrL(y)) = _L e -W)/T + ^__l e -V A (r,T)/T _ (3) 



The singlet and adjoint channels into which it decomposes follow after projection on 
the corresponding representation matrices to be jSj 

e- v ^T)/r = _L (Tr Lt(x)L(y)>, 

e- yA(r ' T)/T = ^^(TrLt( x) TrL(y)) - ^ _ - (Tr ti(x)L(y)). (4) 

The limit Nt — > oo takes us to zero temperature, at which V av ,Vi,VA represent the 
ground state potentials. At finite N t viz. T, they are instead Boltzmann weighted 
sums over the excitation spectrum and interpreted as free energies. Of the above 
expressions only the average correlator is gauge invariant. In |H] it was shown that one 
can complete the singlet correlator gauge invariantly by dressing the Polyakov lines 
with some functional of spatial links, fi[£/j], which under local gauge transformations 
transforms as 

n 9 (x) = g(x)£l(x)tf(t). (5) 
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Here h(t) is an undetermined SU(N) matrix which may be different in every timeslice. 
Replacing L(x) in the above expressions by 

L(x) = fi f (x, 0)L(x)ft(x, 0), (6) 

renders the singlet and adjoint correlators gauge invariant as well. Of course, Q(x) may 
be interpreted as a gauge fixing function which is local in time (such as e.g. Coulomb 
gauge), in which case the dressed correlators are equivalent to the original correlators 
in a fixed gauge. We stress however that this does not affect the exponential decay. 
As was shown in [UJ |U] , the Wilson loop and the correlator of gauge fixed temporal 
Wilson lines decay with the spectrum of the same transfer matrix, only the matrix 
elements are different. One can use this to re-express Eqs. (j3J) by the manifestly gauge 
invariant periodic Wilson loop. Choosing the gauge fixing functions in (Tr L^(x)L(y)) 
to represent the spatial Wilson line (or "string") between x and y, f2(x, t)f2t(y, t) = 
£/(x, y; t), one obtains the periodic Wilson loop, 

(TrLt(x)L(y)) = (Tr Lt(x)C/(x, y; 0)L(y)C/t(x, y; N t )) = W{\x - y\, N t ), (7) 

which is equivalent to the singlet correlator Eq. (jlj) in axial gauge, t/(x, y) = 1. This 
operator is cheap to compute, manifestly gauge invariant and the one we shall use in 
our numerical investigation. However, our observations are valid for any other choice 
of Q(x) which is local in time and transforms as Eq. (j5J), as well. 

The original "average", "singlet" and "adjoint/octet" labeling of the correlators 
refers to the transformation properties of the operators when explicit fields for the 
static sources are introduced. At a given time, interpolating operators for a static 
meson in a color singlet and adjoint state are, e.g., 

0(x, y) = ^(x)f/(x, y)^(y), O a (x, y) = ^(x)f/(x, x )T a f/(x , y)^(y), (8) 

with group generators T a . Here xo represents the center of mass coordinate of the 
meson, at which these operators transform as a singlet and adjoint under gauge trans- 
formations, respectively. Integrating over the static fields and using the completeness 
relation for T a T a , one finds for the corresponding correlation functions 

(O(x,y;0)Ot(x,y;iV t )) 
(O a (x,y;0)O a t( x ,y;iV t )> 



i.e. Eqs. (JH) in axial gauge 



oc (TrLt(x)C/(x,y;0)L(y)C/t(x,y;iV t )) 
1 



'OC 



N 2 -l 
1 



N(N 2 - 1) 



(TrLt(x)TrL(y)) (9) 
(Tr L+(x)£/(x, y; 0)L(y)C/t( x , y; N t )) 
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Expressing the Euclidean expectation values as Hamiltonian traces over complete 
sets of states by means of the transfer matrix formalism ^01 [H] , we find 



(O(x,y;0)O t (x,y;iV 4 )) cx Z^fi 
(O a (x,y;0)O a t(x,y;iV i )) cx Z^fr 

where £7 a (x, y) = U (x, x )T a £7 (x , y) , and the transfer matrix is specified in Eqs. ()13|14|) 
below. Note that after integrating out the source fields, the operators in Eq. (jHJ) reduce 
to their pure gauge parts, 

O a/3 (x,y) = ^(x,y), 0^(x,y) = L/^(x,y). (11) 

In addition to transforming at x under the local gauge group, they also transform like 
a fundamental (F)/antifundamental (F) representation at x, y, respectively. When 
acting on the gauge-invariant vacuum, they hence generate states living in the F(x) cS> 
F(y) <g> l(x ) = F(x) <g> F(y) and F(x) <g> F(y) eg) A(x ) sectors of the full Hilbert space, 
respectively. Keeping the sources saturating the indices at x, y in mind, we shall refer 
to these as the sectors of singlet and adjoint states. 

Evidently, the exponential decay of both correlators in Eqs. (fTU|) is governed by the 
spectrum of the transfer matrix, 

(. - ■) oc £ ce-*^, (12) 

n 

and the energy levels are independent of the detailed form of the operators U(x, y) (or 
other gauge fixing functions fi(x)) [E]. Based on this it was concluded that separate 
potentials in the adjoint /singlet channels can be obtained, which then combine to the 
average potential as in Eqs. (|3I4|) . 

However, this conclusion does not hold. As we shall now show, the so-called adjoint 
and average correlators are mislabeled: they receive contributions exclusively from 
singlet states. Even though the operators in the adjoint correlator transform as an 
adjoint at x , the eigenstates contributing to it do not. The crucial observation is that, 
after integrating out the source fields, the transfer matrix in Eq. (JTUJl has acquired 
indices needed for the time evolution of the string operators. This reflects the fact that 
on all timeslices the Gauss law with static sources is imposed rather than the Gauss 
law without sources as in the case of correlators of particle states. It turns out that, 
for both correlators, the states propagated by the transfer matrix are in the singlet 
sector F(x) <g> F(y). These states combine with the operators to give gauge- invariant, 
non-trivial matrix elements. 



^ /3 (x,y)(f iVt ) /37fo ^(x,y) 



U 



x,yj 



(10) 
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3 Projection on the Hilbert space sectors 



This is best seen by using projection operators on the Hilbert space. The transfer 
matrix in Eq. (fTU|) can be decomposed as 

where T is the familiar transfer matrix in temporal gauge acting on the Hilbert space 
in the presence of external charges [TO] ITT]. i- e - To = exp(— aH ) with H the Kogut- 
Susskind Hamiltonian [12] . The operator 

PX = J Dg gi^)gMR[g], (14) 

where R[g] imposes a gauge transformation g(x) on wave functions, projects on the 
sector with a static quark-antiquark pair. More specifically, it annihilates all wave 
functions not transforming as ipp^JJ 9 } = g{pi) p 1 ip 1 s\U]g\^{y) and maps to ip ua (see 
Appendix El). Note that this is the transformation behavior of our singlet (at x ) wave 
functions. On the other hand, the operator projecting on adjoint states is given by 

PXTat = jDg gUx)g,MDi(g(^))R[g], (15) 

where D^ b (g) are the representation matrices of the adjoint representation, 

Dt b (g) = 2Tr^T a gT b ). (16) 

It is well known [T3] that the "average" correlator can be expressed as a quantum 
mechanical trace over a complete set of eigenstates of T , 

^(TrLt(x)TrL(y)) = JL £ tr [f* P^%] = ^tr [T * P F ^] 

a/3 

= ^£lK>^>| 2 e- E " /T =^£e- E " /T . (") 

na/3 n 

where we have used that p F ®^ = N 2 P^^ is a projector. The presence of this projector 
enforces that only singlet eigenstates |n Qj g) transforming like ipa/3 above contribute to 
the correlator. Hence all energy levels contributing to the thermodynamic sum are 
energies of singlet states, with no matrix elements besides the constant 1/N 2 . To 
understand how this fits together with the supposed decomposition into singlet and 
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adjoint contributions, we similarly rewrite Eqs. (jl(Jj) . Using periodicity of the trace, 
one finds 

(O(x,y;0)Ot(x,y;JV t )> cc Z^fr [f N * Mx,y)L^(x,y) Pjg ] 

= jp ElKI^(^y)^(x,y)|n, Q )| 2 e~ E - /T 

n 

(O a (x, y; 0)O1»(x, y; iV t )) oc Z" 1 tr [if ^ 5 (x, y)^( X) y ) 

= jj* E l(^ 7 |C> 7 a ,(x,y)C/:;(x,y)|^ a )| 2 e"^. 

(18) 

These expressions reveal two non-trivial features: firstly, even in the "adjoint" cor- 
relator projection is onto the F(x) ® F(y) ® l(x ) sector of Hilbert space, while the 
adjoint projector, Eq. (fT5|). does not appear at all. Hence, the energy levels _E n con- 
tributing to both expressions are the singlet ground state potential and its excitations. 
Second, the indices of the operators are contracted with those of the eigenfunctions of 
the Hamiltonian, rather than among each other. This means that there are non-trivial 
matrix elements depending on the operators U (x, y) (or other gauge fixing functions 
fi(x)). As we shall see in Sec. these matrix elements are responsible for the structure 
observed in previous simulations in the literature E] • 



4 Spatial exponential decay 

The correlators Eqs. (|3I4|) can of course equally be viewed as correlation functions in 
space, and often one is interested in their spatial decay with the separation r = |x — y| 
between the static sources. In this case one performs an analogous analysis by defining 
a spatial transfer matrix, and corresponding Hamiltonian, which propagate states along 
the axis defined by the charges. In this case one finds for the average correlator the 
well known result 

^(Tr L+(x)Tr L(y)> = Z' X T r [T L ~ r Tr i,t( x± ) f r Tr L(y ± )\ 

n 

where xi stands for the coordinates perpendicular to the axis x — y. 

The transfer matrix T = TqP acts on the sector of gauge invariant eigenstates of the 
spatial Hamiltonian. The eigenvalues E n depend on the finite N t . In particular, for 
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large N t the theory is confining, and E n ~ N t are energies of torelonic states, while for 
small N t they are gauge invariant screening masses [T^j . 

For the singlet correlator some care has to be taken, because in general a gauge 
fixing function Q(x) is non-local in the correlation direction and thus prohibits the 
definition of a positive transfer matrix. However, our choice of axial gauge is simply 
the spatial equivalent of temporal gauge, for which a positive transfer matrix exists 
and the analysis can be performed. We find 

^(Tr L+(x)l/(x, y; 0)L(y)^( X; y; Nt)) = z^Tr ^L^^L^)] (20) 

^ ^El^l^(x±)|0)| 2 e-^. 

n 

In the second equality, we have used the transformation behavior (j2~B*j) of U, applied 
to the Polyakov line L, which is now transverse to the direction in which the transfer 
matrix acts. Note that the energy eigenvalues in this case are those of the transfer 
matrix T , describing the sector with charges propagating spatially. Hence, with respect 
to the spatial decay there is a difference between the two channels of correlators. This 
difference is also reflected in perturbation theory, where to leading order Eq. (jl9j) is 
dominated by two gluon exchange, compared to one gluon exchange for Eq. (J20|) HE] ■ 



5 Numerical results 

In this section we present the results of our numerical study of pure SU(2) gauge theory 
in 2+1 dimensions. Simulations in that theory are cheap and large Polyakov loops can 
be computed without recourse to error-reducing techniques. In order to unambiguously 
extract the lowest energy eigenvalues and their corresponding matrix elements, the 
number of timeslices Nt has to be large enough for the ground state to dominate the 
exponential decay (i.e. to approach zero temperature, T = l/(aN t ) — > 0). We work 
at gauge coupling (3 = 9, which is large enough for the physical spectrum to be close 
to the continuum, and spatial volumes V = 24 2 known to be free of finite size effects 
[mini- By considering lattices with different temporal extent N t = 24, 26, 28, 30, 32, we 
are able to monitor whether the correlators are indeed described by single exponential 
decay, and cleanly separate exponents and matrix elements. 

The result of this procedure for r/a — 1, 4 is shown in Fig. Q where all three corre- 
lation functions are plotted against N t . Table [T] shows the resulting slope parameters, 
the static potentials E\{r), for r/a = 1, 2, 3, 4. One finds that the potential is the same 
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Figure 1: The three different correlation functions Eqs. (E^p, constructed from dressed 
Polyakov lines, forr = a (left) andr = 4a (right); the lines represent single exponential 
fits. 



r/a 


"average" 


X 2 /dof 


singlet 


xVdof 


"adjoint" 


x7 dof 


WL 


1 


0.09731(2) 


0.61 


0.09749(1) 


0.26 


0.0925(51) 


0.94 


0.0976(3) 


2 


0.1518(7) 


0.51 


0.1526(2) 


0.42 


0.144(6) 


0.56 


0.1529(3) 


3 


0.1923(18) 


0.43 


0.1934(11) 


1.07 


0.1865(72) 


0.32 


0.1935(4) 


4 


0.2312(83) 


1.13 


0.2286(33) 


1.42 


0.239(25) 


0.98 


0.2280(4) 



Table 1: Fitted static potentials aEi(r) from "average", singlet and "adjoint" correla- 
tion functions, as well as the Wilson loop. The values agree within errors. 

for all three, and that it corresponds precisely to the zero temperature singlet potential 
one obtains from a Wilson loop calculation at the same lattice spacing [Bj. 

The matrix elements of the ground states, cf. Eq. ()12)) . are shown in Fig. |2l While 
our signal gets quickly noisy with growing distance, one can clearly discern that in the 
"average" channel there is only the trivial constant 1/N 2 , in accord with our spectral 
decomposition. In the singlet channel on the other hand, there is a non-trivial matrix 
element deviating significantly from one with growing separation, and similarly for the 
"adjoint" channel. 

In the light of this, how can the previous numerical results in the literature be 
understood? If we just extract the potentials according to the prescription of formulae 
Eqs. (I3I4|) . we obtain the plot shown in Fig. El This indeed reproduces previous results 
which were suggestive of three different potentials. However, from our analysis it is 
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Figure 2: T7ie matrix elements of the ground state in all correlators, Eq. Mfy). 



quite clear that this structure is due to the matrix elements, which get exponentiated 
when formulae Eqs. ()3I4|) are used. In order to verify this claim, let us truncate the 
spectral decompositions of the potentials to the ground state (since we are near the 
infinite N t , or T = limit), 



-V av (r)/T 
-Vi(r)/T 
-V A (r)/T 



cfe _i?l(r)/T , 



c^r) e 



Ci r e 



-£i(r)/T 
~Ei(r)/T 



(21) 

(22) 
(23) 



with the same energy E\{r) for all channels, as predicted. Expression (|17|) for the 
decomposition of the "average" potential and relation (j3J) then imply 



1 



cf(r) 



cj(r) 



TV 2 



(24) 



Thus, the would-be adjoint and average channel potentials are really given in terms of 
the singlet static potential Ei(r) and its matrix element 



(25) 



We have inserted the central values of our data for Ei(r), c\(r) in these formulae, Fig- El 
which indeed reproduce the curves extracted by means of the old definition. 
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Figure 3: The potentials computed according to Eq. Q). The "adjoint channel" rep- 
resents a combination of the singlet potential and its matrix element according to 
Eqs. \2'J\) and \2$, which is shown by the full circles. 



6 Conclusions 

By using projection operators on Hilbert space sectors as well as a numerical analysis of 
the zero temperature limit, we have shown that Polyakov loop correlators in each of the 
"average" , singlet and "adjoint" channels defined in the literature receive contributions 
exclusively from singlet states. Their exponential decay in the temporal direction is 
thus determined by the singlet potential and its excitations only. As expected, in 
the zero temperature limit the ground state singlet potential is identical with the one 
extracted from Wilson loops. Correspondingly, at finite temperature one may extract 
only a singlet free energy. At fixed lattice spacing, its correct definition as a Boltzmann 
weighted sum of exponentials without matrix elements follows from Eq. (|17|) to be 

e -fi(r)/T = ^ e -iVT = (TrLt(x)TrL(y)). (26) 

n 

Any other definition, as well as the other channels considered in the literature, expo- 
nentiate (operator dependent) matrix elements, thus faking an r and/or T-dependence 
which is not shared by the physical states. 

Finally, let us note that we do not dispute the existence of an adjoint sector, 
F(x)(g>F(y)(g>A(xo), in the Hilbert space. However, it is not probed by any combination 
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of Polyakov loops, but requires a different operator, which must not close through the 
time boundary in order to truly couple to adjoint states. For a non-perturbative inves- 
tigation, the only way of saturating an open adjoint index is then by a corresponding 
adjoint source. However, in this case the resulting object is a gluelump formed by an 
adjoint source coupled to two fundamental ones in the adjoint representation, which is 
no longer the situation of interest (and considered by perturbation theory), with only 
two fundamental sources in the adjoint. It appears that more work is needed to resolve 
these questions. 

Note added. In a somewhat different use of language, it has been suggested to in- 
terpret hybrid potentials (i.e. angular momentum excited states with two fundamental 
sources as octet potentials, since they appear to approach the gluelump spectrum 
in the short distance limit [18j. According to their transformation behavior, all these 
belong to the singlet sector, since the total system of static quark, antiquark and glue 
is singlet under gauge transformations. While relevant for hybrid mesons, these poten- 
tials are therefore not useful for describing intermediate color-charged states as they 
appear in some nonrelativistic QCD computations These latter states have to 

carry an extra adjoint charge. We thank G. Bali for correspondence on this point. 



A Projectors on charge sectors 

The transfer matrix T acts on the full Hilbert space of all square integrable wave 
functions, including ones transforming nontrivially under gauge transformations. This 
Hilbert space splits into orthogonal sectors with charges in arbitrary representations at 
any lattice point. These are characterized by their transformation properties under the 
local gauge group and can be obtained by acting with appropriate projection operators. 
The projector on gauge-invariant states (with no charges) can be written as 

P = j DgR[g], (27) 

where R[g] performs a gauge transformation with gauge function g, R[g]ip[U] = ip[U 9 ]. 
Specifically, a string between fundamental charges at x, y transforms as 

R\g)U a p(x,y)Bj\g] = ^(x)t/ 7 ,(x,y)^(y) . (28) 

The projector on states with a fundamental charge at x and an anti-fundamental one 
at y is thus given by 

pF^F = N 2 f ^Tr^(x)Tr^(y)%] . (29) 
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The operator P^®^ defined in Eq. (jT3j) maps the component ipp^ of the representation 
F ® F to the component ip au and annihilates all other components and charge sectors. 
The components are here defined by the transformation property 

R{9] *l>Pix = Vv W x )&Uy) • ( 30 ) 

A direct computation shows 

^^ = ^<WWVw (3i) 

and that P^Z, annihilates all other charge sectors: 

Pit (1 - P m ") = • (32) 

ri o o r- 



The normalization of PEgZ, has conveniently been chosen such that 



J2 P X^ = ^u- (33) 

The projection operators N 2 P^®^ (no sum) provide a decomposition of p F ® F y 

P m * = N*Y,PlZ- (34) 

a/i 
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